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. The phase diagram of strongly interacting matter is explored as a function of temper- 
ature and baryon number density. We investigate the possible simultaneous formation 
of condensates in the conventional quark-anti-quark channel (breaking chiral symmetry) 
and in a quark-quark channel leading to color superconductivity: the spontaneous break- 
ing of color symmetry via the formation of quark Cooper pairs. We point out that for 
two massless quark flavors a tricritical point in the phase diagram separates a chiral sym- 
metry restoring first order transition at high densities from the second order transition 
at high temperatures. Away from the chiral limit this tricritical point becomes a second 
order phase transition with Ising model exponents, suggesting that a long correlation 

a.: 



length may develop in heavy ion collisions in which the phase transition is traversed at 



^ ■ the appropriate density. 



^ ■ 1. Towards the phase diagram 

sh ; 

The behavior of QCD at high temperature and baryon density is of fundamental interest 
and has applications in cosmology, in the astrophysics of neutron stars and in the physics 
of heavy ion collisions. Over the past years, considerable progress has been achieved in our 
understanding of high temperature QCD, where simulations on the lattice and universal- 
ity arguments played an essential role. Recently, nonperturbative renormalization group 
methods have been established which account for both the low temperature chiral per- 
turbation theory domain of validity and the domain of validity of universality associated 
with critical phenomena JTJj. They may help to shed some light on the remaining pressing 
questions at finite temperature, like the nature of the phase transition for realistic values 
of the quark masses. 

Our knowledge of the high density properties of strongly interacting matter is rudimen- 
tary so far. There are severe problems to use standard simulation algorithms at nonzero 
chemical potential on the lattice due to a complex fermion determinant. Different nonper- 
turbative methods like the Exact Renormalization GroupQ or Schwinger-Dyson equations 
seem to present promising alternatives. However, finding suitable nonperturbative 
approximation schemes often relies on some knowledge about propagators or the relevant 



1 First investigations in this direction done by D. Jungnickel, C. Wetterich and myself (|] have been 
presented by D. Jungnickel, these proceedings. 
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Figure 1. Phase diagram as a function of T and /x in GeV. The solid curves are second 
order phase transitions; the dashed curve describes the first order transition. 



degrees of freedom. As a first step it seems well justified to consider models which al- 
low us to describe likely patterns of symmetry breaking and to make rough quantitative 
estimates. In this talk, I will present an exploration of the phase diagram for strongly 
interacting matter, using a class of models for two flavor QCD in which the interaction 
between quarks is modelled by that induced by instantons. The model is discussed using a 
mean field approximation. The results are based on recent work done in collaboration 
with Krishna Rajagopal. 

QCD at high density is expected to have a rich phase structure. In addition to the 
nuclear and quark matter phases a number of possibilities like the formation of meson 
condensates or strange quark matter have been discussed. A particularly interesting 
possibility is the formation of Cooper pairs of quarks, where an arbitrarily weak attractive 
interaction between quarks renders the quark Fermi surface unstable and leads to the 
formation of a condensate. Since pairs of quarks cannot be color singlets, a diquark 
condensate breaks color symmetry. For recent reviews of this rapidly growing field of 
research see Pfl. Different attractive channels may lead to a (simultaneous) formation 
of different condensates. In particular, the vacuum of QCD already has a condensate 
of a quark-anti-quark pair. An important ingredient for the understanding of the high 
density phase structure is the notion of competing condensates. One expects, and we 
find ||, that the breaking of color symmetry due to a (ipip) condensate is suppressed by 
the presence of a chiral condensate. Likewise, we find chiral symmetry restoration to be 
induced at lower densities by the presence of a color superconductor condensate. This 
behavior can be quantitatively understood from the effective potential which I discuss 
below. The competition between condensates may greatly simplify calculations. We find 
that previous treatments J7|,[§, which discuss the chiral and superconducting condensates 
separately and inspired this work, are a good approximation once the phase boundary is 
known. It is encouraging that this finding has also recently been confirmed in calculations 
in an instanton liquid model |§. 

The phase diagram which we uncover has striking qualitative features, several of which 
we expect to generalize beyond the model which we consider. The phase diagram is shown 
in Fig. p] for zero quark mass as a function of temperature T and chemical potential \i 
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for (net) quark number. At low temperatures, we observe Q chiral symmetry restoration 
via a first order transition between a phase with low baryon density and a high density 
phase with a condensate of quark-quark Cooper pairs in color antitriplet, Lorentz scalar, 
isospin singlet states.^ There are coexisting (ipip) an d (V^) condensates in this phase in 
the presence of a current quark mass. We find color superconductivity at temperatures 
T < T C A where T C A is of order tens to almost one hundred MeV. The transition we find 
at T C A is second order for the present model, but the order of this transition may change 
once gauge field fluctuations are taken into account. 

In the chiral limit, we find a second order chiral transition at zero chemical potential, 
and a treatment of the present model which went beyond mean field theory, e.g. along 
the lines presented in [ED], would find critical exponents characteristic of the three dimen- 
sional 0(4) universality class. The point in the phase diagram where the second order 
transition meets the low temperature first order transition locates the tricritical poinQ 
The temperature is estimated in our model to be T tc ~ 101 MeV and a chemical poten- 
tial /i tc — 211 MeV. This point has been observed in different phenomenological models 
[p^| . ^5| ,[4|,p!6| . We note that the critical behavior in the vicinity of this point is governed 
by universality (pointed out independently in [[HJ). If two-flavor QCD has a second order 
transition at high temperatures and a first order transition at high densities, then it will 
have a tricritical point in the same universality class as our model. The physics near the 
tricritical point is described by a 6 field theory and there are three independent critical 
exponents which are given quantitatively by our mean field analysis and by logarithmic 
corrections to scaling [T^J. A nonzero quark mass m explicitly breaks chiral symmetry 



and the second order phase transition above T tc turns into a smooth crossover. A small 
quark mass cannot eliminate the first order transition below T tc . Therefore, whereas we 
previously had a line of first order transitions and a line of second order transitions meet- 
ing at a tricritical point, with m^Owe now have a line of first order transitions ending at 
an ordinary critical point. The situation is precisely analogous to critical opalescence in a 
liquid gas system. At this critical point, one degree of freedom (that associated with the 
a) becomes massless, while the pion degrees of freedom are massive since chiral symmetry 
is explicitly broken. Therefore, this transition is in the same universality class as the three 
dimensional Ising model. 

From many studies of QCD at nonzero temperature, we are familiar with the possibility 
of a second order transition, with infinite correlation lengths, in an unphysical world in 
which there are two massless quarks. It is exciting to realize that if the finite density 
transition is first order at zero temperature (see also ||) then there is a tricritical point in 
the chiral limit which becomes an Ising second order phase transition in a world with chiral 
symmetry explicitly broken. In a sufficiently energetic heavy ion collision, one may create 
conditions in approximate local thermal equilibrium in the phase in which spontaneous 
chiral symmetry breaking is lost. Depending on the initial density and temperature, when 
this plasma expands and cools it will traverse the phase transition at different points in the 
(/i, T) plane. Our results suggest that in heavy ion collisions in which the chiral symmetry 



2 There are also indications of a color 6, Lorentz axial vector, isospin singlet condensate which is many 
orders of magnitude smaller than the condensates we treat. 

3 Note that the presence of a tricritical point is precisely what is expected to occur at /j = at a particular 
value of the strange quark mass |[l0|-[l2| . The two tricritical points are continuously connected H . 
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breaking transition is traversed at baryon densities which are not too high and not too 
low, a very long correlation length in the a channel may be manifest even though the 
pion is massive @,|16|- Recently, a number of distinctive signatures have been proposed 



[ [T3|j that could allow an identification of the predicted critical point. 

2. Effective potential 

Our results are obtained from a class of models for QCD where the fermions interact 
via the instanton induced interactions between light quarks. The interaction reflects the 
chiral symmetry of QCD: axial baryon number is broken, while chiral SU(2)l x SU(2)r 
is respected. Color SU(3) is realized as a global symmetry. We note that with the help 
of appropriate Fierz transformations the instanton interaction can be decomposed into 
two parts, where one part contains only color singlet fermion bilinears and the other part 
contains only color 3 bilinears: Si = Si + Sj 3 ^ with (see [|J) 

+0(rf) [ip, V>; -p]Otf) [V>, V>; p] + O a {ao) [ip, -p\O iao)a [ip, p] } 

sf c) = g 2 tJ2 / 7 ^{-oS^;p]0( s)q [V';p] + oJ5^;^0( P ).[V';p]} 

where we have generalized the interaction to allow G\ and G 2 to take on independent 
values. Here p = (2n7rT,p) due to the bosonic nature of the fermion bilinears. The 
bilinears 0( CT ), 0? n \, and Of ao ^ with a = 1, 2, 3 carry the quantum numbers associated 
with the scalar isosinglet (<r), the pseudo-scalar isotriplet (n), the pseudo-scalar isosinglet 
(r/) and the scalar isotriplet (ao), respectively. Similarly, the bilinear (O^), with 
the color index a = 1, 2 or 3, carries the quantum numbers of the color antitriplet scalar 
(pseudo-scalar) diquark. For the a and for the scalar diquark they read 

nez J \ Z7T ) 

where C denotes the charge conjugation matrix and we have supplemented the bilinears 
by suitable form factors F(q) to mimic the effects of asymptotic freedom. The results we 
quote in our exploration of the phase diagram are obtained using the smooth form factor 
F(q) = A 2 /((f 2 +A 2 ) with A = 0.8 GeV and with G\ fixed by requiring a constituent quark 
mass of 400 MeV at/i = T = m = 0in order to obtain a reasonable, albeit qualitative, 
phenomenology and G2 = 3Gi/4 which is motivated in ||. We note that the qualitative 
features which we address do not depend on this specific choice^. 

We note from the signs in Sj that if we choose the sign of G\ such that the interaction 
in the a channel is attractive, so that chiral symmetry breaking is favored, we may expect 



4 There is a phenomenological upper bound which G2 must satisfy. For G2 > 2Gi color would be 
spontaneously broken in the vacuum. 
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Figure 2. The zero temperature thermodynamic potential Q (in GeV 4 ) as a function of 
(j) ~ ^/>Vy at A = (left panel) and as a function of A ~ at = (right panel) for 

several chemical potentials. The curves correspond to (top to bottom) \i = 0, 0.292, 0.35 
GeV. 



condensation in the it and scalar diquark channels also. In the chiral limit, one can always 
make a rotation such that there is no it condensate. A condensate in the pseudoscalar 
diquark channel would break parity spontaneously, but this seems not to be favored by the 
interaction (p. We will use the model to explore condensates in the a and scalar diquark 
channels, i.e. we consider possible simultaneous condensates, (O^) ~ (tpip) for chiral 
symmetry breaking and (0? s \) ~ ( , 0'0) f° r breaking of color symmetry 577(3) — ► SU(2). 

The appropriate tool to study the phase structure of the model is to consider the effec- 
tive potential Q (generating functional of 1PI Green functions at zero momentum). Using 
standard techniques |3| we compute the mean field effective potential ft(4>, A; //, T) as a 
function of two 'order parameters^ and A. Extremizing Q leads to coupled gap equa- 
tions, whose solutions <fto and Ao are related to the chiral condensate and the condensate 
of Cooper pairs by <f) = 2Gi(i)ij^, A = 2G J 2 ('0'0)- At its extrema it corresponds to 
the thermodynamic potential, related to the energy density e, the entropy density s, the 
quark number density n and the pressure P by 

Q((p , A ;/j,,T) = e-Ts-fm = -P . (3) 

As the temperature, the chemical potential or the quark mass changes, the effective 
potential can have several local minima in the ((f), A) plane. Only the global minimum 
corresponds to the lowest free energy state and is favored. The analysis yields the phase 
structure presented in Fig. |l| 

To discuss the competition between condensates let us concentrate on T = 0. Fig. |2| 
shows 'slices' of Q as a function of <fi at A = and as a function of A at = for 
several chemical potentials. We find that f2(0 = 0, A) shows a nonzero minimum Ao 
for arbitrarily small fi: Since the instanton interaction provides an attractive interaction 
in the color 3 scalar diquark channel it renders the Fermi surface unstable and would 

5 We note that the two flavor diquark condensate is not a gauge invariant order parameter. In this 
respect it is analogous to the electroweak sector of the standard model. Likewise, there is the possibility 
of no sharp transition but rather a crossover with qualitatively different physical behavior and the same 
symmetries. 
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Figure 3. The chiral gap <po (solid line) and the diquark gap A (dashed line) as a function 
of the chemical potential \x in GeV for zero quark mass m (left panel) and for m — 10 
MeV (right panel). 

lead to the formation of a diquark condensate. However, one observes that the chiral 
condensate is favored for /i < 0.292 GeV and the diquark condensate vanishes identically. 
For larger \x chiral symmetry restoration occurs via a first order transition, while the 
diquark condensate jumps to a nonzero value in the high density phase. We therefore 
have the situation that the competition with the chiral condensate leads to a first order 
color superconductor transition which would be continuous otherwise. One also notes 
from Figs. § that without considering the possibility of diquark condensation, i.e. for 
an identically vanishing A, chiral symmetry restoration would have occured at larger 
chemical potential when Q(<fi, A = 0) crosses the origin. In the chiral limit we observe a 
strong competition in the sense that where one condensate is nonzero the other vanishes 
(see also Fig. |3|). 

It is instructive to consider the physics away from the chiral limit. A nonzero quark mass 
breaks chiral symmetry explicitly even in the high density phase. Indeed, for an average 
current quark mass of 10 MeV one observes in Fig. ^| simultaneous condensates with 
comparable magnitude immediately upon the completion of the transition. Comparison 
with the chiral limit (left panel) shows that the diquark condensate is not significantly 
disturbed by the presence of <p . This insensitivity can be understood by noting that 
although the Cooper pairs have low momentum, they are formed from quarks which have 
momenta close to the Fermi surface. Adding a quark mass m <C \x does not significantly 
affect the density of states or the interactions of the quasiparticles with momenta of order 
fi. The phenomenon of color superconductivity seems quite robust. Even away from 
the chiral limit it is not significantly disturbed by the simultaneous presence of a chiral 
condensate in the high density phase as long as m + 0o is smaller than the chemical 
potential. 

3. Outlook 

We have presented here a first model calculation of the QCD phase diagram that takes 
into account the strong competition between the chiral and the color antitriplet scalar 
diquark condensate. The latter involves only the light u and d quarks of two colors. 
Though the strange quark is much heavier than the two light quarks, it may not be 
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neglected for chemical potentials much larger than m s or typical scales of a few hundred 



MeV. There is a compelling symmetry breaking scheme [M, involving three massless 
flavors, which locks color and flavor to a residual global SU(3) symmetry. It is not clear 
what phase will be realized with physical mass it, d and s quarks. As pointed out in ||, 
one may think of having the two flavor condensate first and color-flavor locking at higher 
density. A first investigation may be performed along the lines presented here [[| but 
finally more sophisticated methods than the mean field analysis will be needed to settle 
these questions. A promising possibility may be the use of truncated nonperturbative 
flow equations which have been successfully applied to finite temperature in the past. 

The phase diagram which we uncover in our two flavor study has striking qualitative 
features, most remarkably the presence of a tricritical point. For nonzero quark masses a 
situation similar to the liquid-gas nuclear transition arises which has been much studied 
in low energy heavy ion collisions. We observe the (tri) critical point to emerge as a result 
of a second order transition/crossover in one region of phase space (high T, small /z) and 
a first order transition for the same order parameter in another region (high /i, low T). 
The question of whether the tricritical point is realized in QCD is closely connected with 
the still unsettled question of the order of the high temperature (/i = 0) transition, and 
has to be addressed in a three flavor study. If the strange quark mass is too small, or if 
the axial U(l) symmetry is effectively restored about the transition, then we may have a 



first order transition at high T which is driven by fluctuations [2C]. In this case, a line of 



first order transitions connects the T and /i axes. However, one should note that if the 



finding [19] of a continuous high density crossover for physical mass u, d and s quarks is 
realized in nature, then again a critical point in the phase diagram may emerge: A high 
temperature first order transition line ending in an Ising endpoint in the high density 
region. The involved long-range correlations in any of these scenarios would be appealing 
both from the theoretical and the experimental point of view. 
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